Weak Lensing of the Primary CMB Bispectrum 
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The cosmic microwave background (CMB) bispectrum is a well-known probe of the non- 
Gaussianity of primordial perturbations. Just as the intervening large-scale structure modifies the 
CMB angular power spectrum through weak gravitational lensing, the CMB primary bispectrum 
generated at the last scattering surface is also modified by lensing. We discuss the lensing modifica- 
tion to the CMB bispectrum and show that lensing leads to an overall decrease in the amplitude of 
the primary bispectrum at multipoles of interest between 100 and 2000 through additional smooth- 
ing introduced by lensing. Since weak lensing is not accounted for in current estimators of the 
primordial non-Gaussianity parameter, the existing measurements of /nl of the local model with 
WMAP out to Zmax ~ 750 is biased low by about 6%. For a high resolution experiment such as 
Planck, the lensing modification to the bispectrum must be properly included when attempting to 
estimate the primordial non-Gaussianity or the bias will be at the level of 30%. For Planck, weak 
lensing increases the minimum detectable value for the non-Gaussianity parameter of the local type 
/nl to 7 from the previous estimate of about 5 without lensing. The minimum detectable value of 
/nl for a cosmic variance limited experiment is also increased from less than 3 to ~ 5. 
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PACS numbers: 98.70.Vc,98.65.Dx,95.85.Sz,98.80.Cq,98.80.Es 



I. INTRODUCTION 

The weak lensing of cosmic microwave background 
(CMB) anisotropy angular power spectrum is now well 
understood in the literature [l|, Q . The modifications re- 
sult in a smoothing of the acoustic peak structure at large 
angular scales and an increase in power below a few ar- 
cminute angular scales corresponding to the the damping 
tail of CMB anisotropics Q . 

The angular power spectrum of the lensing potential 
out to the last scattering surface can be established with 
quadratic estimators that probe lensing non-Gaussianity 
at the 4-point level of a CMB map Q . Such a reconstruc- 
tion of the lensing potential is helpful for CMB B-mode 
studies of polarization especially in the context of 
searching for the signature of the primordial tensor modes 
0. This is due to the fact that in addition to inflation- 
ary gravitational waves, the B-modes of CMB polariza- 
tion also contains a signal generated by lensing of scalar 
E-modes with a peak in power at a few arcminute angu- 
lar scales The lensing reconstruction has now been 
applied to the existing Wilkinson Microwave Anisotropy 
Probe (WMAP) data leading to a ~ 2cr to 3cr detection 
of gravitational lensing in the CMB through a correlation 
between the reconstructed lensing potential and tracers 
of the large-scale structure such as radio galaxies [1, [§] . 

In parallel with the progress on lensing studies with 
the CMB, the search for primordial non-Gaussianity us- 
ing the CMB bispectrum with constraints on the non- 
Gaussianity parameter /nl is now an active topic in cos- 
mology 0, d, [il- The 5-year WMAP data is consis- 
tent with — 9 < /nl < 111 at the 95% confidence level 
for the local model [l3], though a non-zero detection of 
primordial non-Gaussianity at the same 95% confidence 
level with 26.9 < /nl < 146.7 is claimed elsewhere using 
the WMAP 3-year data [ll|- This result, if correct, has 
significant cosmological implications since the expected 



value under standard inflationary models is /nl ^ 1 
[H, [H, [13, [H, [H, [13, mi , though alternative models of 
inflation, such as the ekpyrotic cosmology [H, [131, gener- 
ally predict a large primordial non-Gaussianity with /nl 
at few tens. 

Just as the CMB power spectrum is modified by lensing 
from potential fluctuations of the intervening large-scale 
structure [H, the CMB bispectrum will also be modi- 
fied by gravitational lensing. The correlation between 
the projected lensing potential and CMB secondary ef- 
fects, such as the integrated Sachs- Wolfe (ISW) effect 
or the Sunyaev-Zel'dovich (SZ) effect, leads to a non- 
Gaussian signal at the three point level [1^, [2^ . These 
secondary non-Gaussianities are expected even if the pri- 
mordial perturbations are Gaussian and impact existing 
primordial non-Gaussianity parameter measurements by 
introducing a small, but unavoidable, bias [27l. [28l. [29| . 

Beyond secondary non-Gaussianities, weak lensing by 
the intervening large-scale structure maps the intrinsi- 
cally non-Gaussian CMB sky to a different anisotropy 
pattern when observed today. Thus the bispectrum one 
reconstructs with a CMB map, assuming it to be of the 
expected form at the last-scattering surface due to pri- 
mordial non-Gaussian perturbations, will result in a bi- 
ased estimate of the primordial non-Gaussianity param- 
eter. The existing estimator can be modified to account 
for lensing modifications and to obtain a bias-free esti- 
mate of the non-Gaussianity, but at the expense of fac- 
torizability that has allowed fast computation of the bis- 
pectrum in existing analyses [30| . Since lensing modifies 
the anisotropy pattern by smoothing the fluctuations, a 
change in the minimum detectable non-Gaussianity pa- 
rameter /nl for a given experiment is expected to be 
different from the existing values in the literature [l0| . 

In this paper, we present a general derivation of the 
lensed CMB primary bispectrum and quantify above 
statements on the changes imposed by lensing for the 
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detection of primordial non-Gaussianity. We find that 
the non-Gaussianity parameter measured with WMAP, 
ignoring lensing, will result in an estimate of /nl for the 
local model that is biased low by about 6%, when mea- 
surements are extended out to Imax ^ 750. Furthermore, 
with lensing, the minimum detectable level of /nl with 
Planck is increased by roughly 40% from less than 5 to 
about 7 and the cosmic variance limit of /nl is increased 
from 3 to 5. 

This paper is organized as follows: we first discuss the 
CMB primary bispectrum of the local type in Section II. 
Some basic ingredients related to the lensing calculation 
is presented in Section III. We derive the lensing effect 
on the bispectrum, under both flat-sky and all-sky for- 
mulations, in Section IV. We discuss our results and con- 
clude with a summary in Section V. In illustrating our 
results we make use of the standard flat ACDM cosmo- 
logical model consistent with WMAP with Qb = 0.042, 
fla = 0.238, h = 0.732, = 0.958 and r = 0.089. 



II. 



CMB PRIMARY BISPECTRUM 



The CMB temperature perturbation on the sky, 
Q{n) = Ar(n)/T, is decomposed into its multipole mo- 
ments 



e(n) = ^e,,„rr(fi)- 



(1) 



The angular power spectrum and bispectrum of CMB 
anisotropics are defined in the usual way, respectively, as 



^1 h h 
nil W3 



(SiimiBi^m^efjTOj) - ( ]Bi^i^i^, (2) 



where, for the bispectrum, we have introduced the 
Wigner-3j symbol (see the Appendix of Ref. [l^ for some 
useful properties of this symbol) . 

The CMB bispectrum is generated by a coupling of the 
local-type with a quadratic correction to the Newtonian 
curvature such that 

$(x) = $i(x) + /nl [$i(x) - (<i>i(x))] (3) 

where $ l (x) is the linear and Gaussian perturbation and 
/nl in the non-Gaussianity parameter, which is taken to 
be scale independent p^ . 

In Fourier space, we can decompose equation Q as 



with 



$(k) =$L(k) + /NL$NL(k), 



$NL(k) = I ^$4k + ki)<I>2(ki 

d^ki 



(4) 



where is the linear power spectrum, defined as 

($(k)$(k')) = (27r)3,5(k + k')P4.(fc) . (6) 

The multipole moments of the anisotropy can be writ- 
ten as 



47r(-i)' 



(27r)3 



Hk)gTiimi,{k), (7) 



where $(k) from above is the primordial curvature per- 
turbation in the Fourier space, and gxiik) is the ra- 
diation transfer function. With <I'L(k) and <I>nl, the 
moments can be separated into two components with 

^Im — + ^Im ■ 

The CMB angular power spectrum can be defined us- 
ing the transfer function and the power spectrum of dom- 
inant linear fluctuations as 



(8) 



Q - - / k'dkP^{k)g^i{k). 



Using the definition of the angular bispectrum (equa- 
tion [2|), the primordial temperature anisotropy bispec- 
trum can be written as 
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which can be simplified to 
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BiM, = 201,1,1, / r^drbi^,{r)b^,{r)b^^^{r) (10) 
Jo 

+btir)bfi'^{r)btir) + bl\r)bf,,{r)bf,{r), 

where 

2 r°° 

hf{r) = - / edkP^ik)gTiik)ji{kr), (11) 
2 f°° 

bi'-i.r) = /nl- / k^dkgTiik)jiikr), (12) 

7^ Jo 



and 

Ghhh 



(2;i + l)(2Z2 + l)(2?3 + l) fh h h 
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(13) 

When illustrating our results, we make use of a modi- 
fied code of CMBFAST [HI for the standard flat ACDM 
cosmological model to fully calculate radiation transfer 
functions when generating the CMB primary bispectrum. 



III. WEAK LENSING BASICS 

The effects of weak lensing can be encapsulated, under 
the Born approximation, in the radial projection of the 
gravitational potential ($), given as p32l | 



/f/31, d A(r — r'^ 

j^P.iki), (5) 0(fi) ^ -2 1 ^^'4r^<J>(r(n),r'), (14) 
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FIG. 1: The CMB bispectrum for the equilateral case {h — 
h ~ h ~ I) with (solid line) and without (dashed line) lensing. 
Here we plot Z*_Bf;;/(27r)^ as a function of the multipole I for 
one of the sides. We assume /nl ~ 1. The lensing effect can 
be described as a decrease in the amplitude of the bispectrum 
when I < 1700 with an increase at higher multipoles. 

where r(z) is the hne-of-sight comoving distance (or look- 
back time) to a redshift z from the observer with last 
scattering surface at = r{z = 1100), and dA{r) is the 
comoving angular diameter distance. In a spatially flat 
universe, (Ia — ^ r. Here, we ignore the time-delay effect 
as it is small compared to the geometric lensing effect 
captured by equation pij) [s^. 

The calculation related to the CMB bispectrum de- 
scribed below requires the angular power spectrum of 
lensing potential (f), which can be decomposed into the 
multiple moments as 

0(n)=^0z™l^™(fi), (15) 

Im 

with the lensing power spectrum defined using 

{4>l'm4>l'm') = 5l,l>5m,rn'Cf tO obtain 3 

Cf = ^ j edkP^{k)[ll^''{k)]'' , (16) 

where 

iy'™(r) = -mr) '!iZ~"\ - (17) 
dA[rjdA(rs) 

Here F{r) describes the radial evolution of potential fluc- 
tuations. Modifications to the CMB anisotropics, gener- 
ated at higher order in lensing potential fluctuations, are 
at the level of at most 5% relative to those due to the lens- 
ing potential angular power spectrum [33 |. The bispec- 
trum of lensing potentials, due to the non-linear evolution 
of density perturbations, also modifies the CMB primary 



bispectrum, but these changes can also be ignored since 
the lensing potential bispectrum is at the order (^Cf^ , 

while changes we describe are first order in the angu- 
lar power spectrum of the lensing potential. Using the 
Limber approximation, equation (|16[) can be further sim- 
plified, but we do not make use of this approximation in 
numerical calculations illustrated here since the flat-sky 
form of the potential power spectrum is known to bias 
lensing results of the power spectrum by about 10% at 
all multipoles [1]. 

IV. LENSING OF THE CMB BISPECTRUM 

We first give a treatment of the lensing of the CMB bis- 
pectrum assuming a flat-sky approximation and discuss 
a derivation under the spherical sky later. 

A. Flat-sky Case 

Weak lensing deflects the path of background photons 
resulting in a remapping of the observed anisotropy pat- 
tern on the sky. Following an approach similar to [3[ , we 
write the lensed temperature anisotropy as 

e(fi) = e[n + V9i(n)] 

« e(n)-HV,0(fi)v^e(fi) 

+ iv,(/)(n)V,(/.(n)V'V^e(n) + ... (18) 

Here, 0(n) is the unlensed CMB temperature anisotropy, 
0(n) is the lensed anisotropy, and V0(n) is the lensing 
deflection angle for the CMB photons. 

Taking the Fourier transform, as appropriate for a 
fiat-sky, we write the lensed temperature anisotropy in 
Fourier space as 

e(li) = J dne(n)e-''i-" 

= e(ii)- J ^e{i,')Lih,h'), (19) 

where 

L(li,li') = ,^(li - li') (li - li') • li' 

x(li".li')(li-li'-li")-li' (20) 

to the second order in lensing potential in the perturba- 
tive expansion. 

The observed angular power spectrum of CMB 
anisotropics under weak lensing is discussed in Q. The 
resulting power spectrum consists of both the unlensed 
intensity and a perturbative correction related to the 
lensing effect. Making use of the expansion and after 
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FIG. 2; The squeezed configurations {h ^ h ^ h) of the CMB bispectrum with (soUd hue) and without (dashed fine) lensing. 
The left panel is for h — 10 and right panel is h = 100. We vary I — h with Is = I + h m both cases and plot l'^ Bij^^i+i^/{2tv)'^ 
as a function of I. Again, we take /nl = 1- In these configurations, the lensing effect can be described as an overall decrease 
in the amplitude of the bispectrum when / < 1200. This suggests that lensing by the intervening large-scale structure leads to 
a less non-Gaussianity in the CMB map. 



some straight forward calculations, we obtain the lensed 
anisotropy power spectrum as 



where 



l^R) 



4tt 



dl P Cf 



(21) 



(22) 



Here, R describes the variance of the deflection angle. For 
ACDM cosmology, 0rms = VR ~ 2.6'. This derivation 
makes use of the flat-sky approximation to describe the 
lensing effect on CMB anisotropy power spectrum. When 
the expressions derived in the previous section for Cf and 
cf under the exact spherical-sky treatment are used in 
equation (PT|) . the lensed CMB power spectrum can be 
derived with a less bias than using, say, the flat-sky result 
for cf in the same expression Q. 

Keeping the flat-sky approximation, we can define the 
angular bispectrum as 

(e(ii)e(i2)e(i3)) ^ {2^f5{\, + h + h)Bf,^,^,.^,) , (23) 

and following the approach similar to the lensed angu- 
lar power spectrum that led to equation pT|) . the lensed 
bispectrum can be expressed as 



-^(1142,13 
-"(1142,13 



(24) 



(il + q + il) 

m(12-1') 



R 



-f^(ll,l2-l',l3+l 



+-^(li-lM2+lM3)(^i 



+B' 





(Il-|-l',l2,l3-l 



')(l3 



2 

l'(ll + 12 - 1') • 1' 

1') • l'(l3 + li - 1') • 1' 
1') • l'(l3 + h~ 1') • 1' 



Note that we have identified the flat-sky bispectrum 
as i?^^ to distinguish from the all-sky bispectrum 

(2 h ■ "^^^ related through 
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h h h\ J{2li + l){2l2 + l){2h + l) 
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(ii, 12,13) 
(25) 



B. All-sky Treatment 

The derivation related to the lensing of the CMB bis- 
pectrum under the more appropriate spherical sky can 
be obtained by replacing the Fourier components with 
spherical harmonic multipole moments. In this case, the 
lensed field can be represented as Q 



e 



Irn 



e 



dnY,"*V.0(n)V*e(n) 
dnl^™*V,0(ri)Vj0(n)V''V^e(n) 



0;™ + ^ ^ ei",„" (26) 

I'm' l"m" 

^ Tmm' m" , ^ \ ^ -tmm'm"m"' 

^ll'l" "T ^ 2^ 'Pl"'m""-'ll'l"l"' 



V" m'' 



where, the integrals over the spherical harmonics were 
replaced, in the last step, by the geometrical factors 



(27) 
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Using equation ([2]), the lensed CMB temperature bis- 
pectrum can then be expressed as 



leading to 



7712 1713 



B 



E 



(28) 
(29) 



mim2m3 



^1 h h 
mi 1712 TO3 



(0iimi0i2m2©;3 



"2 E E E i'^hmiQl2rn2Ql'^7n'^ 

"E E E E ^^'i^i^'i'K^'s™: 



h'^m'^<f>i'^m'Jli^V^q +2 Perm. 



Noting that the Wigner-3j symbol obeys the identity 

^1 k h \ f h h I'i \ ^ hiV^^niim'^ 

2-^ I mi 7712 7713 i I ?77i m2 m(j' ) ^ (2I3 + 1) ' 

mim2 / ^ / 

(30) 

we can re-write the lensed bispectrum as 

^z!m3 = <M3 + ^<^2^3E^?^i+2Pe™. 

'3 

+ E^feE^/t^2 + 2Perm. (31) 



'■2 '3 ^—-1 '2 



where 



^ —1— y y j"3™>3m^ 

(2/3 + 1) ^ ^ i3i3i3i3 

"13 

^ ^ f h h h \ f h I2 ^3 

^ 2-^ 2-^ I mi m2 m3 I I mi m2 mg 

m!,mi,'7n'J mim2m3 



^ J^2^2^^2 ^^i3^i2"^3 
^2^2 '2 '3 '2 '3 



The sum of the geometric term j™3™3™3m3 Qygj. 

o I3I3I3I3 J 

yields [1] 



47r 



(33) 



Hence, with summation over m3 leading to a factor (2^3 + 
1), the expression for Si simplifies to 

1 9/' + 1 

^1 = -3^3(^3 + 1)^^(^^ + 1)^^. (34) 



In order to evaluate S2, we first re-express I;™™ ™ 



^ll'l" 



' m m' m" ' ' 



where Q 



z r I" 




r(r+i)+r'(r+i)-/(z+i) 

Then the expression for S2 can be re-written as 
S2 = //2;y^' //aiy.!,' E E 



(35) 



(36) 



h h I3 
mi m2 m3 



h I2 '3 ) ( '2 I2 I2 ) ( '3 ^2 ^3 

mi m2 777.3 / V '^2 7n2 m'2 j \ 7773 m2 m'^ 

fl2l'2l2 /i3i2i3 (^1)' ^'^^'M ^1 r ' (37) 



where, in the last step, we have introduced the Wigner-6j 
symbol . The values of the Wigner-6j symbol can be 
computed numerically with a fast and efficient recursive 
algorithm [11]. 

Finally, substituting the expressions for Si and S2 in 
equation (|3ip and including all permutations in a single 
expression, we can write the lensed bispectrum as 



b: 



^ 1 ^2^3 



1 - ni hih + 1) + hih + 1) + ^3(^3 + 1) 



B 



E^/ 



£ f / 1 \n J ^1 ^2 ^3 
Jl2lplhlq\-^} ] I q p 
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p I q 

h h I3 

q p I 
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where 



21 + 1 

An 



Cf 



(38) 



(32) and n = {I + p + q). 



V. RESULTS & DISCUSSION 

We illustrate the modification to the equilateral con- 
figurations of the bispectrum in Figure 1, where we plot 
l'^ B^i / {2tt)'^ as a function of the multipole /. The primary 
CMB bispectrum assumes /ml = 1 and is calculated with 
the full radiation transfer function gTi{k). The lensing 
description makes use of the all-sky treatment to calcu- 
late both cf and the lensed bispectrum. The flat-sky 
expression gives a result consistent with the all-sky ex- 
pression at better than 5% at all multipolcs if the all-sky 
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FIG. 3: The signal-to-noise ratio for a detection of the CMB bispectrum with (thick lines) and without (thin lines) lensing. 
The long-dashed lines show the case for WMAP and dot-dashed lines for Planck. The left panel shows the signal-to-noise ratio 
as a function of ^3, while the right panel shows the cumulative signal-to-noise ratio below I3 in the x-axis. Note the overall 
reduction in the signal-to-noiso ratio (when ^3 ~ 1500) in the case of lensing relative to the case where lensing is ignored. 



expression for C[ is used in both calculations. The differ- 
ence is at the level of 10% if the two expression make use 
of the two separate calculations of Cf , as in the case of 
the angular power spectrum In the case of the equi- 
lateral configurations of the bispectrum, the lensing effect 
can be best described as a smoothing and a decrease of 
the amplitude of non-Gaussianity power in the equilat- 
eral configuration of the bispectrum when I < 1500 and a 
subsequent increase in the bispectrum amplitude at small 
angular scales. 

In Figure 2, we show two squeezed configurations 
{h ^ I3 ^ h) of the bispectrum, with the short length 
fixed at cither 10 (left panel) or 100 (right panel), as 
a function of the multipole I of one side with the third 
side fixed at either I + 10 or I + 100. Without lensing, 
a comparison of Figure 1 and Figure 2 reveals a well- 
known result in the literature that the local type of the 
CMB bispectrum is dominated by squeezed configura- 
tions with one small side and two large sides for the 
bispectrum mode shape in the multipole space. With 
lensing, the amplitude of the squeezed configurations is 
significantly reduced when two of the sides have lengths 
I > 1200 in the multipole space. This can be again de- 
scribed as a smoothing effect with lensing by the inter- 
vening large-scale structure; removing or "washing out" 
the primordial non-Gaussian signature in the CMB map 
at angular scales below a few arcminutes. Thus, when 
lensed by the large-scale structure, the primordial non- 
Gaussian CMB sky appears more Gaussian at arcminute 
scales when studying the non-Gaussianity at the three- 
point level. At same angular scales, however, the CMB 
sky appears more non-Gaussian due to lensing at the 
four-point level probed by the trispectrum (sBj . 

The removal of the non-Gaussianity is associated with 
the squeezed configurations, which dominate the overall 
signal-to-noisc ratio for the detection of the primary bis- 



pectrum without lensing. Although non-Gaussianity is 
reduced for the squeezed configurations, lensing leads to 
an increase in the amplitude of the bispectrum for equi- 
lateral configurations where ~ ^2 ^ ^3 > 1500. This 
increase, however, is insignificant in terms of the overall 
signal-to-noisc ratio as the contribution to the cumulative 
signal-to-noisc coming from these configurations is lower, 
owing to the higher variances associated with foregrounds 
and instrumental noise at these angular scales. 

To further quantify this statement, we plot, in Fig- 
ure 3, the signal-to-noise ratio calculated as 



E 



Q^tot ^tot ^tot 
^1 ^2 ^3 



(39) 



where the noise variance calculation involves all con- 
tributions to the angular power spectrum with C*°* = 
Ci + Cf^"^ + Ni where we include the lensed CMB power 
spectrum (C/), secondary anisotropics {Cf'^'^), and the 
noise power spectrum (Ni) for both WMAP and Planck. 
For secondaries, we include the SZ power spectrum calcu- 
lated with the halo model [13 and make use of the noise 
calculations from Ref. [2^ for WMAP and Planck. For 
the case involving an experiment limited by the cosmic 
variance, we set Ni = 0. In the left panel of Figure 3, 
we plot d{S/N)^/dl3 as a function of ^3, while in the 
right panel we plot the cumulative signal-to-noise ratio 
out to ^3 in the x-axis. In the case where lensing is not 
included, signal-to-noise ratio estimates for the bispec- 
trum detection are consistent with previous calculations 
in the literature [l0| . 

With lensing, however, the signal-to-noise ratios are 
changed. As can be seen from the left panel of Figure 3, 
there is an overall reduction in the signal-to-noise ratio 
when I3 1500. This difference comes from the pre- 
viously described decrease in the amplitude of the non- 



7 



1000 




1000 



100 



1000 



FIG. 4: Contour plots of d{S/N)^ /dloglma^dloglmin (equa- 
tion |40]) as a function of Zmax and Zmin for the primary bispec- 
trum (top panel) and the lensed primary bispectrum (bottom 
panel). We take /nl = 1. 



Gaussianity in squeezed configurations of the bispectrum 
with lensing imposed. To further understand the differ- 
ences in the signal-to-noise ratio of the lensed primary 
bispectrum, we plot, in Figure 4, the quantity 



d{S/Nf 



d\0g lmaxd\ogl„ 



max*- mill 



E 



l = ln 



g^tot ^tot^tot 



(40) 

with two separate estimates for and to esti- 
mate this quantity without and with lensing, respectively. 
Note that the overall signal-to-noise ratio comes from in- 
tegrating this quantity over the variables ^max and /min 
and we include the factor /max^min to account for the 
logarithmic scaling. A comparison of the two panels in 
figure 4 reveals an overall decrease in the amplitude at 
I ^ 10^ in the case with the lensed primary bispectrum 
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FIG. 5; Contour plot of the difference d[(5/Ai')fcnsGd ~ 

(S'/iV)unlcns;,d)/o'log^maxrflogimin aS a fuUCtloU of ^max and 

Imin (same as the difference between bottom and top panels 
of Figure 4). 
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FIG. 6: The bias in the primordial non-Gaussianity parameter 
/nl when lensing is ignored in the estimation. (-!-)/(-) shows 
the sign of A/nl//nl between the true and measured values 
(see text for details). 



relative to the primary bispectrum alone. 

In figure 5, we plot the difference of the two as a con- 
tour plot to show that lensing results in an overall de- 
crease in the signal-to-noise ratio in the squeezed configu- 
rations when /,nax 10^ and laiin < 10^, while there is an 
increase in the signal-to-noise ratio when Imax ^ 3 x 10"^ 
for all values of Imin- These plots demonstrates the same 
trends described with respect to figure 2 involving a de- 
crease in the amplitude of the squeezed configurations of 
the lensed bispectrum. While there is a slight increase 
in the lensed bispectrum amplitude at I3 > 2000, such 
small angular scales are not probed by Planck. Even in 
the cosmic variance limit, unfortunately, this small in- 
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crease is insignificant given that at these same angular 
scales secondary anisotropics dominate the bispectrum 
noise variance. In terms of the cumulative signal-to-noisc 
ratio values shown in figure 3, the minimum /nl to detect 
the bispectrum with Planck and a cosmic variance lim- 
ited experiment is increased by about 30% to 40% from 
/nl 5 for Planck to /nl ^ 7. The cosmic variance- 
limited detection threshold for /nl is increased from 3 to 
5. 

While the difference in cumulative signal-to-noise ratio 
seems insignificant for an experiment like WMAP, weak 
lensing could impact existin g rn easurements of the non- 
Gaussianity parameter /nl jisl . [T^ . To understand the 
lensing bias introduced to /nl, we follow the discussion 
in Ref. [1^ and note that the current estimators of the 

non-Gaussianity parameter use /nl = Uli , 

with 

pq 

where B^^^^i^ is the primary bispectrum. and Cpq is the 
covariance matrix for bispectrum measurements involv- 
ing triplets of p = (hhh) and q = (^'i^2^3). This esti- 
mator is the optimal estimator for non-Gaussianity mea- 
surements, but given the complications associated with 
estimating the covariance, existing studies make use of 
a sub-optimal estimator which approximates the covari- 
ance with variance C~^, w cr~'^{li,l2,h)Saa' , and intro- 
duces a linear term to equation (j4ip to minimize the vari- 
ance of /nl [13 ■ Note that N is the overall normalization 
factor that can be calculated from equation ([1T|) by re- 
placing B°^^ with B^. 

While weak lensing modifies the observed bispectrum 
^f^hlh ^ /nl^z^j^/j, existing measurements make the 
assumption that Bf^^^^i^ = JnlB^i^i^- This results in 
a biased estimate of /nl from the true value of the 
non-Gaussianity parameter /nl°- The fractional differ- 
ence of this bias A///nl = (/nl° " /nl)//nl can be 
calculated through the covariance between the lensed 
and unlensed CMB primary bispectrum A///nl = 

1 - E/?z!z.Z3'^;'A!z.Z3/I:5z?z.Z3^"'^;!m3]' ^l^cre we 
have simply written the variance as cr ^. We plot 
A///nl as a function of / to which non-Gaussianity pa- 
rameter measurements are performed in figure 6. Ex- 
isting measurements with WMAP data probe out to 
^max 750 and wc find that existing estimates of /nl 



arc biased by ~ 6%. For Planck, if lensing is ignored, the 
bias is at the level of 30%. 

This bias is not the same fractional difference in the 
signal-to-noise ratio that one can infer from figure 3 
since the fractional difference in the signal-to-noisc ra- 
tio with and without lensing involves a ratio of the form 
[EBj%,a-^Bf^i,,Jj:Bf,i,,,a-^Bf^^,^^J. Note that in 
future an unbiased estimate of /nl can be obtained by 
replacing -Bj^j^is equation (|4T|) with the lensed bispec- 
trum B^^i^i^ and recalculating the normalization factor 
N with the lensed primary bispectrum. Unfortunately, 
while without lensing the CMB primary bispectrum of 
the local model factorizes into two separate integrals with 
bf" and bf^ (described in Section II) , this factorizabil- 
ity is no longer preserved when lensed and impacts an 
easy estimation of the non-Gaussianity parameter with 
the existing estimator [s^l. For Planck and other CMB 
experiments that can probe down to small angular scales 
for primordial non-Gaussianity measurements, it will be 
necessary to implement an estimator that accounts for 
the lensing effect. 

To summarize the main results of this paper, we have 
discussed the primary CMB bispectrum generated at the 
last scattering surface, but observed today after it is weak 
lensed by the intervening large-scale structure. Unfortu- 
nately, as we have found, weak lensing leads to an overall 
decrease in the amplitude of non-Gaussianity with the 
biggest change on the squeezed configurations of the bis- 
pectrum that dominate the overall signal-to-noise ratio 
when studying the primordial non-Gaussianity parame- 
ter. For an experiment such as the Wilkinson Microwave 
Anisotropy Probe (WMAP), the modifications imposed 
by lensing results in an estimate of /nl of the local model 
that is biased low by about 6%. For a high resolution ex- 
periment such as Planck, the lensing modification to the 
bispectrum must be accounted for when attempting to 
estimate the primordial non-Gaussianity. The minimum 
detectable value of /nl for a cosmic variance limited ex- 
periment is ~ 5. 
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